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Abstract - Recent structural findings have shown that dynamin, a cytosol protein playing a 
key-role in clathrin-mediated endocytosis, inserts partly within the lipid bilayer and tends to self- 
assemble around lipid tubules. Taking into account these observations, we make the hypothesis that 
individual membrane inserted dynamins imprint a local cylindrical curvature to the membrane. This 
imprint may give rise to long-range mechanical forces mediated by the elasticity of the membrane. 
Calculating the resulting many-body interaction between a collection of inserted dynamins and a 
membrane bud, we find a regime in which the dynamins are elastically recruited by the bud to form 
a collar around its neck, which is reminiscent of the actual process preempting vesicle scission. This 
physical mechanism might therefore be implied in the recruitment of dynamins by clathrin coats. 

endocytosis / clathrin / dynamin / membrane inclusions interactions 

Resume - Des donnes structurales recentes out montre que la dynamine, une proteine du cytosol 
qui joue un role cle dans I'endocytose clathrine-dependante, s'insere partiellement dans la bicouche 
membranaire et tend a s'auto-assembler autour de tubules lipidiques. En tenant compte de ces 
observations, nous faisons I'hypothese que les dynamines impriment localement une courbure cylin- 
drique dans la membrane. Cette empreinte pent engendrer des forces elastiques de longue portee. 
En calculant I'interaction multi-corps entre un ensemble de dynamines inserees dans la membrane 
et une capsule endocytotique, nous trouvons un regime dans lequel les dynamines sont recrutees 
elastiquement par la capsule pour former un collier autour de son con, ce qui rappelle le processus 
precedant la scission des vesicules d'endocytose. Ce mecanisme physique pourrait done etre implique 
dans le recrutement des dynamines par les capsules de clathrine. 

endocytose / clathrine / dynamine / interactions entre inclusions membranaires 



I. INTRODUCTION 

In eukaryotic cells, membranes of different organelles are functionally connected to each other via vesicular 
traiisport. Formation of transport vesicles from invaginated buds of the plasma membrane is called endocyto- 
sis [ij . In clathrin-mediated vesiculationjVesicle formation starts with the assembly on the donor membrane 
of a highly organized "coat" of clathrins [J , which acts both to shape the membrane into a bud and to select 
cargo proteins 0, 0, 0, 0, • The mechanism by which an invaginated clathrin-coated bud is converted to a 
vesicle (scission) involves the action of a cytoplasmic GTPase protein called dynamin A ^l. D ynamins form 
oligomeric rings at the neck of deeply invaginated membrane buds and induce scission |9l llfl||. How exactly 
dynamin is recruited and how the scission actually occurs remains unclear |llL Il2j . In this paper we propose 
that dynamin recruitment by clathrin coats could be driven by long-ranged physical forces mediated by the 
membrane curvature elasticity. 

Crvp-electron microscopy has recently revealed the detailed structure of the clathrin coats at 21 A resolu- 
tion Clathrin units, also called "triskelions" , have a star-like structure with three legs. Initially solubilized 
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FIG. 1: Schematic representation of a cytosol protein partly inserting within a lipid bilayer and inducing a local 
membrane curvature via a binding region. 



into the cytoplasmic fluid, they self-assemble onto the membrane surface into a curved, two-dimensional solid 
scaffold. The latter is a honeycomb made of hexagons and pentagons (geometrically providing the curvature) 
the sides of which are built by the overlapping legs of the clathrin triskelions. In the plasma membrane, 
clathrins usually interact with "adaptor" transmembrane proteins, which also serve to select cargo proteins. 
However it has been shown that clathrin coats can readily self-assemble onto protein- free liposomes 0, Q] . 

Dynamin is known to be solubilized in the cytosol as tetramers, and to aggregate in low-salt buffers into 
rings and spirals Dynamin also self-assembles onto lipid bilayers, forming helically striated tubules that 
resemble the necks of invaginated buds (tube diameter c± 50 nm) 10] . Addition of GTP induces morphological 
changes: either the tubules constrict and break jl^, or the dynamin spiral elongates 0|. These findings 
suggest that the scission of clathrin-coated buds is produced by a mechanochemical action 0, . 

At the earlier stages of the budding process, dynamins already strongly interact with bilayer membranes. 
Indeed, in vivo studies showed that dynamin binds acidic phos p holi pids in a way that is essential to its 
ability to form oligomeric rings on invaginated buds 0, 0, I21L l22|. Using a model lipid monolayer 
spread at the air-water interface, it was shown that dynamins actually penetrate within the acyl region of 
the membrane lipids |2,l| . This finding was recently confirmed by the three-dimensional reconstruction of 
the dynamin structure by cryo-electron microscopy at 20 A resolution [23 | : dynamins form T-shaped dimers 
the "leg" of which inserts partly into the outer lipid leaflet. 

It was long ago suggested |2^ I23 that particles inserted within bilayers should feel long-range interactions 
mediated by the elasticity of the membrane. Indeed, a protein penetrating within a bilayer and binding its 
lipids — such as dynamin — may in general produce a local membrane curvature (see Fig. ^ . Because of the 
very nature of the curvature elasticity of fluid membranes, this deformation relaxes quite slowly away from 
its source, and the presence of another inserted particle produces an interference implying an interaction 
energy [23. This holds as long as the separation between the inclusions is smaller than the characteristic 
length ^o- = k/ct, where k ~ 60 k^T is the bending rigidity of the membrane and a is the tension of 
the membrane. At separations larger than ^o-, the membranes flattens out and the interaction vanishes 
exponentially. Note that the membrane tension is not a material constant like k; it is an effective force per 
unit area, which is most probably biologically regulated and is usually of the order of 10^^ to 10~^ times 
the surface tension of ordinary liquids l28j. At this point, anticipating on our model for the clathrin- 
dynamin system, let us state that we shall formally assume u = in this paper, which amounts to assuming 
that the relevant distances between the inclusions (i.e., the distance between the neck of the clathrin bud 
and the dynamins) are smaller than ^g.- This means that our model should rather apply to weakly tense 
membranes, e.g., tr ~ 10~^mJ/m^, for which ^o- ~ 400 nm (which is quite larger than the typical size of 
the clathrin buds ~ 80 nm). Although the actual value of a in the vicinity of clathrin buds is unknown, 
such a small tension agrees with recent measurements on biological membranes (erythrocytes membranes 
interacting with their cytoskeleton) In the case of stronger tensions, we expect our results to hold 

nonetheless, the dynamins being "captured" when their Brownian diffusion brings them at a distance from 
the bud less than ^o-- 

The first detailed calculation of the membrane mediated interaction was performed for two isotropic 
particles each locally inducing a spherical curvature pol l30l | . The interaction was found to be repulsive, 
proportional to the rigidity k, of the membrane and to the sum of the squares of the imposed curvatures; it 
decays as i?^^, where R is the distance between the particles. The case of anisotropic particles producing 
non-spherical membrane deformations is even more interes ting, sin ce their collective action on the membrane 
is expected to have nontrivial morphological consequences l3lll32Ll33l |. The local deformation of a membrane 
actually involves two distinct curvatures, associated with two orthogonal directions (as in a saddle or in a 
cylinder). Recent calculations showed that the interaction between two anisotropic inclusions is very long- 
ranged and decays as j^ls^ls^. It is always attractive at large separations and favors the orientation 
of the axis of minor curvature along the line joining the particles |35|- Note that these elastic interactions 
prevail at large separations, since they are of much longer range than other forces, such as van der Waals or 
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FIG. 2: Piece of a model membrane showing a clathrin coated bud and the imprints of inserted dynamins. 
screened electrostatic interactions. 



II. MODEL 

Among the above informations, let us outline the three points that are essential for our model, (i) Clathrin 
coats are solid scaffolds that rigidly shape extended parts of the membrane into spherical caps, (ii) Dynamins 
are solubilized proteins that partly insert within the membrane bilayer. (iii) Inserted membrane hosts that 
imprint a local membrane curvature interact with long-range forces of elastic origin. 

We therefore build the following model. We consider a clathrin coated bud as being a membrane patch 
bearing a constant, fixed spherical curvature. Technically, we shall build the bud by placing a large number 
of point-like spherical curvature sources at the vertex of a hexagonal lattice (see Fig. [SJ . Within the present 
formalism, this is the simplest way to define a rigid, almost non-deformable, spherically curved zone. Since 
the dynamins will not penetrate the bud, we expect that our results will not depend on whether the bud 
is geometrically enforced (which would be conceptually simpler but technically harder here) or built by an 
inclusion array. Note that the shape of the neck around the clathrin bud will result from the minimization 
of the total energy, and will therefore not be enforced artificially. 

Because dynamins partly insert within the membrane and seem to accommodate cylindrical curvature, 
we model them as sources locally imprinting a cylindrical curvature. We place a large number of such 
"dynamins" on a membrane in the presence of an artificial bud as described above (Fig. and we study 
whether the latter will recruit or not the dynamins through elastic long-range forces. 

A. Long-range elastic interactions between many membrane inclusions 

The elastic interaction between N isotropic or anisotropic membrane hosts can be calculated from first 
principles I^tI l3^ Is^ . The membrane is described as a surface which is weakly deformed with respect 
to a reference plane. Without this assumption, analytical calculations are virtually impossible. An obvious 
consequence is that we can accurately describe only weakly invaginated buds; nevertheless, we expect that 
our results will hold qualitatively for strongly invaginated buds. The membrane hosts are described as 
point-like sources bearing two curvatures ci and C2, associated with two orthogonal directions. These values 
represent the two principal curvatures that the hosts imprint on the membrane through their binding with 
the membrane lipids (assuming the binding region is itself curved). For instance, a spherical impression 
corresponds to ci/c2 = 1, a cylindrical impression corresponds to Ci/c2 = 0, and a saddle- like impression 
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corresponds to Ci/c2 = —1. This over-simplified model actually contains the essential ingredients responsible 
for the long-range elastic interactions between membrane inclusions: for protein hosts of a size comparable 
to the thickness of the membrane, it yields accurate interactions for separations as small as about three 
times the particles size. Note that the curvature actually impressed by a particle could be affected by the 
vicinity of another inclusion, we shall however neglect this effect for the sake of simplicity (strong binding 
hypothesis) . 

The point-like curvature sources describing the membrane hosts diffuse and rotate within the fluid mem- 
brane, because of the forces and torques exerted by the other membrane hosts and of the thermal agitation 
fceT. We parameterize the orientation of a particle by the angle 9 that its axis of minor curvature, i.e., the 
axis associated with min(|ci|, |c2|), makes with the x-axis, in projection on the {x,y) reference plane. Given 
N inclusions with specified positions Xi and y^, orientations 9i, and curvatures cu and C2i, for i — 1 ... TV, 
we calculate the shape of the membrane satisfying the N imposed point-like curvatures and we deter- 
mine the total elastic energy of the system. We thereby deduce the TV-body interaction between the hosts 
-Fint(. ■ ■ ,Xi,yi,di,cii,C2i, ■ ■ ■)■ The mathematical details of this procedure are sketched in Appendix IXI 

B. Pairwise interactions 

Before studying the collective interaction between model dynamins and clathrin coats, let us describe 
how point-like spherical and cylindrical sources interact pairwise (in the absence of membrane tension, as 
discussed in Sec. Pi. 

The membrane distortion produced by two inclusions modeled as point-like spherical curvature sources is 
shown in Fig. U^. Each inclusion appears as a small spherical cap away from which the membrane relaxes 
to a flat shape. As evidenced by the plot of the interaction energy (see Fig. such spherical inclusions 
repel one another. Calling k the bending rigidity of the membrane, a the thickness of the membrane (which 
is comparable to the size of the inclusions), c the curvature set by the inclusions and R their separation, our 
calculation gives (see Appendix 

i^el-cl(i?) ^8^KaV (-^)^ (1) 

for the the asymptotic interaction at large separations, in agreement with previous works pfll ls^ . l35l | . Note 
that the plot given in Fig. corresponds to the exact interaction within our model, not to the asymptotic 
expression (^. 

The membrane distortion produced by two inclusions modeled as point-like cylindrical curvature sources 
is shown in Fig.lSb- Each inclusion appears as a small cylindrical cap away from which the membrane relaxes 
to a flat shape. The interaction between two such hosts depends on their relative orientation. The minimum 
energy is found when the axes of the cylinders are parallel to the line joining the inclusions. As evidenced 
by the plot of the interaction energy (see Fig.lSb), the interaction is then attractive. It therefore turns out 
that two such hosts produce a weaker membrane deformation when they are close to one another than when 
they are far apart. As described in Ref. |35j |. when their curvature is strong enough, such inclusions tend to 
aggregate and to form linear oligomers. Their asymptotic interaction energy is given by 

Fdy_dy(i?)^-87rA^aV(|)'. (2) 

It decays as R'^, hence it is of longer range than 

Finally, we show in Fig. |3t the membrane distortion produced by the interaction between a spherical 
source and a cylindrical one. The latter is oriented in the direction that minimizes the energy. As evidenced 
by the plot of Fig. the interaction is attractive at large separations and repulsive at short separations, 
with a stable minimum configuration at a finite distance. Calling c the curvature set by the cylindrical 
inclusion and c' the one set by the spherical inclusion, our calculations give the asymptotic interaction 

Fci_dy(i?) ^ -^nna^cc' (-^)' . (3) 

We may therefore expect that dynamins will be attracted by clathrins coats; however, owing to the non- 
pairwise character of the interaction |35l |. it is necessary to actually perform the corresponding many-body 
calculation. It is also necessary to check whether thermal agitation will or will not disorder the inclusions. 
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FIG. 3: Shape of a membrane distorted by two inclusions imprinting local curvatures and interaction energy as a 
function of separation. Distances are rescaled by the membrane thickness a and energies by na^c? . (a) spherical 
inclusions of curvature c; (b) cylindrical inclusions of curvature c; (c) spherical inclusion of curvature c and cylindrical 
inclusion of curvature 0.2 c. The shapes are calculated from Eq. ljA13fl . the interaction energies from Eq. l|Il2^ . 

III. COLLECTIVE INTERACTIONS BETWEEN MODEL DYNAMINS AND CLATHRIN BUDS 

As described in Sec. El we build a model clathrin-coated bud by placing in a membrane TVci point-like 
spherical inclusions of curvature Cci on a hexagonal array with lattice constant h. Here, we have chosen 
A^ci — 37 and b — 3a. By changing the curvature Cci, we can adjust the overall curvature of the clathrin 
scaffold, thereby simulating the growth of a vesicular bud. We then add A'dy = 40 point-like cylindrical 
sources of curvature Cdy modeling inserted dynamins. 

To study the collective behavior of this system under the action of the multibody elastic interactions (see 
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FIG. 4: Phase diagram representing the typical equilibrium configurations of a system of model dynamins in the 
vicinity of a curved clathrin scaffold. Cdy and Cd are the curvatures associated with the dynamins and clathrins, 
respectively, in units of the inverse membrane thickness a. In region G, the dynamins form a "gas" non interacting 
with the clathrin scaffold. In region L, the dynamins form a system of linear oligomers non interacting with the 
clathrin scaffold. In region R, the dynamins form a ring around the clathrin scaffold, which is reminiscent of real 
endocytosis. 



Sec. Ill A|) and of thermal agitation, we perform a Monte Carlo simulation. The details of the simulations 
are given in Appendix^ To prevent unphysical divergences of the elastic interaction energy, it is necessary 
to introduce a hard-core steric repulsion preventing two inclusions to approach closer than a distance d. 
Since the size of the inclusions imprints is of the order of the membrane thickness a, we have chosen d = 2a. 
Actually, at such microscopic separations, other short-ranged interactions intervene, the details of which are 
still unknown. Here, we disregard them, since our interest lies in the mechanism by which the recruitment 
process and the formation of dynamin collars is driven. In a later stage, which we do not model here, 
dynamin rings are further stabilized by bio-chemical interactions 

The results of the Monte Carlo simulations are summarized in the phase diagram of Fig.^ in terms of the 
curvatures Cd and Cdy of the clathrins and dynamins imprints, respectively. Here, we have chosen to span Cci 
between to 0.2 a~^: for a lattice constant b = 3a and assuming a 40A, this corresponds for the clathrin- 
coated bud to a maximum curvature of radius p ~ b/{a Cd) ~ 60 nm. Since our clathrin patch has 7 spherical 
sources on its diameter, the size of the bud is 7h ~ 85 nm. These values are typical for clathrin-mediated 
endocytosis @. For the dynamins, we have spanned Cdy between 0.1 and 0.4a~^, which corresponds 
to a maximum curvature of the imprint ~ 0.1 nm^^. As for the membrane bending rigidity, we have taken 
K = 60A;bT, since for biological membranes at room temperature n lies between 50 and 100 /cbT" |4(l.l4ll|. 

The phase diagram displays three regimes (see Fig. 0J: a state in which the dynamins are disordered in 
a gas-like fashion (G), a state in which the dynamins form linear oligomers that do not interact with the 
clathrin bud (_L), and a state in which the dynamins form a ring around the clathrin bud (R). In region (i?), 
due to the shallowness of the dynamin imprints, the system is disordered by thermal agitation. Increasing 
the curvature of the dynamin imprints increases the elastic attraction between the dynamins (see Fig Ot) 
and leads to the formation of linear oligomers (L) . These oligomers wrap around the clathrin bud when the 
latter is sufficiently developed (R). Typical snapshots corresponding to the four points a, b, c, d in Fig. 0] 
are shown in Fig. [3 Note that in Fig.Eb the dynamin collar is rather "gaseous" due to the weakness of the 
dynamins' imprints, while in Fig.|3i the ring is tight and well-ordered. 

It is difficult within the present paper to discuss the nature of the boundaries between the different 
"phases" displayed in Fig. 2| Because we are considering a rather small (thermodynaniically speaking) 
number of dynamins, and because they are under the influence of a localized curvature field, the transition 
lines between the various regions of Fig. 21 are actually broad. Attempting to describe them as first- or 
second-order transition lines should not be meaningfuU. 
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(c) 



(d) 



FIG. 5; Typical snapshots showing the equilibrium arrangement of model dynamins (bars) in the vicinity of the 
clathrin scaffold (hexagonal array). The figures (a), (b), (c) and (d) refer to the corresponding points in the phase 
diagram of Fig. 2] 




FIG. 6: Membrane shape corresponding to point (d) in Figs.|l]andEl showing the self-assembly of a ring of dynamins 
around a clathrin bud. 



IV. CONCLUSION 



In this paper, we have shown that if membrane inserted dynamins produce cylindrical imprints and if the 
latter are sufficiently curved, then the resulting long-range forces mediated by the membrane elasticity are 
strong enough to overcome Brownian motion and bring them into a collar around the neck of a clathrin bud. 
Of course, simple diffusion could also bring dynamins around clathrin buds, and their binding into a ring 
could be the result of specific bio-chemical interactions. However, if a cylindrical imprint can speed up this 
process, then evolution may have selected it. 

To test this model, one might look experimentally for linear oligomers of dynamins (see Fig.lSJ^). However, 
since "gaseous" rings are also possible (see Fig. |S|d), the existence of such linear aggregates may not be 
necessary. It would be more interesting to directly check, e.g., by cryo-electron microscopy, the shape of the 
dynamin region that penetrates within the membrane. 

Finally, note that our model is obviously over-simplified: i) many other integral proteins float around 
dynamins, ii) dynamins may interact with various lipidic domains within the bilayer, iii) the membrane may 
have a spontaneous curvature due to its asymmetry, iv) fluctuations are not only thermal but also active, 
and hence could be larger than we estimate, and v) large values of the membrane tension could shorten the 
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range at which the dynamins are recruited (see Sec. P). Nonetheless, we beheve that our model correctly 
captures the effects of the anisotropic elastic interactions. 
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APPENDIX A: MANY-BODY INTERACTIONS BETWEEN POINT-LIKE CURVATURE 

SOURCES 



Let us outline the derivation of the interaction between N anisotropic point-like sources that locally imprint 
a curvature on the membrane. As explained in the text, such constraints modelize a wide class of membrane 
inclusions, including transmembrane and cytosol proteins partly inserted within the membrane. 

For small deformations u{x, y) with respect to the (x, y) plane, the free energy associated with the curvature 
elasticity of a membrane is given by [i^l : 



Fci^'- I dxdy{V^u) 



(Al) 



Indeed, for small deformations, the Laplacian V^u(r) is equal to the sum of the membrane's principal 
curvatures at point r = {x,y). The material parameter k is the bending rigidity. 

Determining the shape of the membrane in the presence of inclusions at positions imprinting local 
curvatures requires minimizing the elastic energy IjAip with local constraints on the membrane curvature 
tensor. In the small deformation limit, the elements of the latter are given by the second spatial derivatives 
of the membrane shape: u^xx{i^), u,xy{T'') and u^yy{r). Introducing 3A^ Lagrange multipHers Afj to enforce 
the curvature constraints, the Euler-Lagrange equation corresponding to the constrained minimization is 



N 

E 



+Ky S,xy (r - r„) + A;;^ S^yy (r - r„) 



(A2) 



where S{r) is the two-dimensional Dirac's delta and a comma indicates derivation. By linearity, the solution 
of this equation is 



3N 

u(r) = ^A^r^(r), 



(A3) 



where the A^'s and F^'s are the 3A^ components of the column matrices 



Al 

yy 

A2 

XX 

V J 



and 



G(r) 



r(r) 



1 

16^' 





n) 


G^xyi^ ^ 


n) 


G,vv(y - 


n) 




r2) 



(A4) 



' hii 



(A5) 



is the Green function of the operator V^V^, satisfying the equation V^V^G(r) = 
We introduce a column matrix K containing the values of the 3A^ constraints 



5{v) 



K = 



/W,x2;(ri)\ 
U^xyi^l) 

U,xx{'r2) 



(A6) 
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With u{r) given by Eq. I|A3(I . enforcing the constraints yields the following equation for the Lagrange 
multipliers: 



3N 



(A7) 



where the 3iV x ?>N matrix M is given by 



M 



mil mi2 
m2i m22 



\mNi 



miN 



(A8) 



in which the rriQ^'s are matrices of size 3x3 defined by 

(G^xxxx (^/3a) G_xxxy {^!3a) G ^xxyy (^/3a)^ 

G^xxxy G ^xxyy (^/3a) G ^xyyy 

G^xxyy {r/3a) G ^xyyy Pa) G ^yyyy (f^a)/' 



where vpa = fa — r^. Setting 



Ya-Yp^ rpa [cOS eia/J X + siu O^p f] , 



(A9) 



(AlO) 



yields explicitely 



^ / cos(4 6'a/3) - 2cos(2 6'q/3) sin (2 6*0/3) [2 cos (2 ^q/j) - 1] -cos (4 6*0/3) \ 

ma/3 = T — — sin (2 6*0/3) [2 cos (2 6*0^) - 1] - cos (4 6*0/3) - sin (4 6^^) - sin (2 6*0-/3) • 

^'^'"/5" V -cos (4 6*0/3) - sin (4 6*0/3) - sin (2 6*0/3) cos (4 ^0/3) + 2 cos (2 6*0/3) / 

(All) 

Integrating Eq. IjAljl by parts and taking into account the constraints yields the elastic energy 

Fei = iKK*M-iK, (A12) 
where K* is the transpose of K. From Eqs. \K?\\ and \kl\ . the equilibrium shape of the membrane is given 

by 

u(r) = K*M-ir(r). (A13) 

When a = /3, mo/3 as given by Eq. I|A11(I diverges: indeed Eq. I|A1|I correctly describes the membrane 
elastic energy only for distances r ^ rg, where is of the order of the membrane thickness. It is therefore 
necessary to introduce a high wavevector cutoff r'^^ in the theory. From the definition of the Green function 
G(r), we deduce, in Fourier space 

G„(r).yA^^. (A14) 
Hence, introducing the cutoff, we obtain 

G^xxxxm = 7^ r cos^ ^ = (A15) 

^0 



(2^)2 7o 32^r2' 



and similarly for the other elements of the matrix l)A9|l . With the above prescription, we obtain 

1 /3 1\ 
^^Trrg 3/ 
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As an illustration, let us consider the case of two identical isotropic inclusions, each prescribing the 
curvature c. Then 

K* = (c,0,c,c,0,c), (A17) 
and, from Eqs. HA12I) . HASp . (|A11|) . and (jA16p . with ri2 = R, the interaction energy is 

512 7rAt(roc)' 
^ci - , . 4 T^Tl ' 



in which we have discarded a constant term. Setting tq — a/2, we indeed obtain the leading asymptotic 
interaction (Q. This special choice for tq allows to match the result of Goulian et al. (1993), which was 
obtained from multipolar expansions. It should be noted that the interaction given by Eq. IjAlSp is exact 
within the present formalism (for r larger than ~ a), whereas multipolar expansions can only give in analytical 
form the leading asymptotic orders. 

When many inclusions are present, the matrix M and its inverse, which determines the interaction energy 
through Eq. HA12|) . can be easily calculated numerically once the positions of the inclusions are defined. 

APPENDIX B: MONTE CARLO SIMULATIONS 

The Monte Carlo simulation that we perform employs the standard Metropolis algorithm For given 
positions and orientations of the particles representing the dynamins, the energy is numerically calculated 
from Eq. (|A12I) . At each Monte Carlo step, we perform a Metropolis move consisting in either a translation 
or a rotation of one arbitrarily chosen dynamin particle. The amplitude of the moves is adjusted in order 
to keep an average acceptance rate of 50%. We confine the dynamins inside a circular box (of radius 80 a) 
centered around the clathrin lattice, which is kept fixed. To take into account the hard-core repulsion (see 
Sec. IIII|) . we simply reject any move bringing two particles closer than the minimum approach distance d 
(here 2a). Note that in this simulation the membrane is not discretized: the interaction energy that we use 
fully takes into account the elasticity of the membrane. 
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